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The polarization of quarkonium states produced in hadron 
collisions exhibits strong non-perturbative effects - for exam- 
ple, at small transverse momentum pt charmonia appear un- 
polarized, in sharp contradiction to the predictions of pertur- 
bation theory. The quark-gluon plasma is expected to screen 
away the non-perturbative physics; therefore those quarko- 
nia which escape from the plasma should possess polariza- 
tion as predicted by perturbative QCD. We estimate the ex- 
pected J /ip polarization at small p t , and find that it translates 
into the asymmetry of the e + e~(/i + /i~) angular distribution 
W{0) ~ 1 + a cos 2 e, with a ~ 0.35 0.4. 



The possibility to form quark-gluon plasma in heavy 
ion collisions is an intriguing problem of strong inter- 
action physics. To establish the formation of plasma, a 
number of signatures were proposed; here we will concen- 
trate on heavy quarkonia. Suppression of heavy quarko- 
nium states has been suggested long time ago by Matsui 
and Satz [1] as a signature of the deconfinement phase 
transition in heavy ion collisions. Their, by now well- 
known, idea is that the Debye screening of the gluon ex- 
changes will make the binding of heavy quarks into the 
bound states impossible or unlikely once a sufficiently 
high temperature is reached. The lack of quarkonium 
states would thus signal deconfinement; this effect was 
indeed observed and studied in detail at CERN SPS by 
the NA38 [2] and NA50 Collaborations [3]. The results 
on J/tp production at RHIC have recently been presented 
by the PHENIX Collaboration [4]. The observations of 
quarkonium suppression have been interpreted as a signal 
of quark-gluon plasma formation [5]. However, different 
conclusions were reached in [6] , where it was argued that 
the effect may arise due to quarkonium collisions with the 
comoving hadrons. Additional tests of the quark-gluon 
plasma formation could help to clarify the situation. 

In this note we propose to use for the diagnostics of the 
quark-gluon plasma those heavy quarkonia which escape 
from it. This would require experimental measurements 
of quarkonium polarization, which can be reconstructed 
from the angular distributions of quarkonium decays - 



dileptons and/or photons. For J/tp states, one would 
need to measure the angular distribution of electrons (or 
muons) in the J/tp — » e + e~ decay in J/tp rest frame 
relative to the direction of its momentum. (We will con- 
centrate on J/tp's at relatively small p tl which dominate 
the total production cross section.) 

Let us first formulate what we mean by the quark- 
gluon plasma, since different definitions sometimes may 
result in misunderstanding. We define the quark-gluon 
plasma as a gas of quarks and gluons in which the in- 
teractions can be described by perturbative QCD and 
non-perturbative effects are either absent or can be ne- 
glected. We will not need to specify the properties of this 
state of matter in more detail to develop our idea; let us 
now turn to the dynamics of quarkonium production. 

It is well-known that the description of the data on 
heavy quarkonium production within the framework of 
perturbative QCD (pQCD) meets with siginificant diffi- 
culties. Both the absolute values of the measured pro- 
duction cross sections of hidden heavy flavor states and 
the relative abundances of different quarkonia are not de- 
scribed well within the perturbative framework, but per- 
haps the most spectacular failure of pQCD is the polar- 
ization of the produced quarkonia. Even an extension of 
a perturbative approach based on non-relativistic QCD 
[7], which allows certain non-perturbative physics, does 
not allow to explain the polarization measurements [8] . 

Meanwhile, the description of heavy flavor production 
in perturbative framework has been largely successful 
(even though there are some problems there as well). The 
reason for this is easy to understand - the production of 
heavy flavors occurs at short time scale ~ l/2mQ, where 
m,Q is the heavy quark mass, whereas the binding of the 
produced heavy quarks into quarkonium is a softer pro- 
cess characterized by the time scale of Ttdnd ~ l/ e > where 
e is the typical binding energy; for a Coulomb interac- 
tion, e ~ OL s rnQ <C rriQ. The binding process is thus far 
more likely to be affected by non-perturbative phenom- 
ena, which manifest themselves both in the magnitude 
of the production cross section and in the polarization of 
the produced quarkonia. 
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Consider now the production of quarkonium states in 
relativistic heavy ion collisions. The typical time scale 
for the production of semi-hard partons with transverse 
momentum k t is r ~ 1/fct; for example, in the gluon satu- 
ration scenario T proc i ~ 1/Qs, where Q s is the saturation 
scale which at RHIC energies is about Q s ~ 1-^2 GeV [9]. 
It is thus likely that while these produced partons will not 
significantly affect the production of heavy quarks (which 
happens at earlier time), they will influence the binding 
of heavy quarks in quarkonia since T proc i < Tund ■ 

High energy density of the produced partonic state 
is expected to result in the destruction of the non- 
perturbative vacuum structure. Indeed, lattice QCD cal- 
culations show that quark and gluon condensate "evap- 
orate" above the deconfinement phase transition [10]. It 
may be expected that non-pcrturbative vacuum fluctua- 
tions are suppressed even if the thermalization does not 
take place - a specific example is given by the suppres- 
sion of instantons in the saturated gluon environment 
[11]. As a result, the processes in this high-density par- 
tonic state of matter should be described by the weak 
coupling, perturbative methods. As a matter of fact, 
as we assumed above, one may define the quark-gluon 
plasma as a collective state of quarks and gluons the 
dynamics of which is governed by perturbative interac- 
tions. Therefore, the formation of heavy quarkonium 
states should also be adequately described by perturba- 
tion theory, and the predictions of pQCD for the polar- 
ization of heavy quarkonia should be vindicated. Dense 
parton matter may then screen out of existence a large 
part of quarkonia, as proposed originally [1], but those 
of them that survive will carry the information about 
the mechanism of their formation throughout the colli- 
sion. Of course, the interactions of quarkonia at the later 
stages of the heavy ion collision may wash out their polar- 
ization somewhat, but relatively small interaction cross 
sections and the heavy quark symmetry, suppressing the 
spin flips of heavy quarks, should prevent quarkonia from 
"forgetting" their initial polarization entirely. 

Let us illustrate this idea in more detail using the ex- 
ample of J/ip polarization. There are two mechanisms 
of J/ip production in hadron collisions - direct, when 
J/ip is produced by perturbative and non-perturbative 
interactions of gluons and quarks, and cascade, when 
J/ip is created as a result of decays of C-even cc states, 
Xc — ► J/ip + 7. In quark-gluon plasma, the cascade pro- 
duction mechanism should be at least as important as 
direct production. Indeed, in the lowest order of pertur- 
bation theory, J/ip is produced by the three gluon fusion 
or by two gluon fusion followed by the gluon emission 
off the cc system. In both cases the probability of J/ip 
production is proportional to oP s {m c ). The probability 
of Xc production is proportional to a s {m c ), i.e. it is of 
lower order in a s , which however is largely compensated 
by the branching ratio B(x2 — ► J/ip + 7) — 20% for the 



J/ip production. 

In hadron collisions the direct mechanism comprises 
typically about 60% of the observed J ftp's (for a review 
of the data, see [12]), which seems to suggest that an es- 
sential part of J/ip production in hadron collisions is of 
non-perturbative origin. Direct calculations confirm this 
conclusion. In rcf [13] J/ip production cross section in 
ttN interactions was calculated in perturbation theory: 
two gluon fusion into cc with the subsequent gluon emis- 
sion (the so-called color-singlet model [14]). The result 
is about 8 times smaller than the data. Similar situation 
holds also for xi production: the calculated cross section 
is by factor of two smaller than the experimental one 
(see [13] for details). Additional mechanism of % 2 pro- 
duction [15] in the framework of the color-octet model [7] 
involves the formation of the color octet cc state which 
then decays by color El transition to X2 ■ Evidently, this 
mechanism perturbatively is suppressed by extra power 
of a s and is essential only if it is non-perturbative. The 
cross section of xi production is very small in perturba- 
tion theory, but noticeable experimentally (xo does not 
contribute substantially to the J/ip production because 
of a small branching ratio of xo —> J/ip + 7 decay - about 
1%). (The contributions from various sources to the J/ip 
production in tt~ N collisions at the incident energy of 
185 and 300 GeV and the results of theoretical calcula- 
tions can be found in [13]; the comparison shows that the 
production of charmonium states in hadronic collisions is 
in an essential way non-perturbative). 

Let us now turn to J/ip polarization as reconstructed 
from the angular distributions of electrons (muons) from 
the J/ip — > e+e~(/i+,u~) decays. Generally the electron 
(muon) distribution has the form 

W{6) ~ 1 + a cos 2 e, (1) 

where is the emission angle of e + (or /i + ) relative to the 
direction of J/ip motion in its rest frame; at small pt, this 
direction coincides with the direction of the beam. The 
value a — 1 corresponds to the transverse polarization, 
a = —1 - to the longitudinal polarization, and a = 
to unpolarizcd J/ip. In perturbation theory, in the case 
when J/ip is produced through the xi ~^ J/^P + 7 decay, 
the coefficient a in Eq. (1) is determined unambiguously 
(at small pt): a = 1 [16]. This comes from the fact 
that X2 is produced by two-gluon fusion, gg — ► xii f° r 
which the effective interaction is f^Q^v, where is 
the energy-mometum tensor of the gluon field and f^ u is 
the wave function of xi ■ Since M „ has only J z = ±2 spin 
projections on the direction of gluon momenta (indeed, 
M „ may be considered as a source of the graviton field), 
the same spin projections has the xi- As a result, J/ip 
produced via X2 decay is transversely polarized, J z = ±1 
and thus a = 1. 

This conclusion is somewhat modified when the initial 
transverse momenta of the gluons are taken into account. 
This reduces the value of a to [16] 
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1 + a 6 2 /r 
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where #g ~ 4(p|)/M^. The average transverse momen- 
tum of gluons is expected to increase with energy and the 
atomic number of the colliding systems. For example, in 
the gluon saturation scenario pt ~ Q s ~ A 1 ^s x / 2 , with 
A ~ 0.25; at RHIC energies in central Au — Au collisions 
Q s ~ 1 GcV [9]. For p t ~ 1 GcV, the formula Eq.(2) 
yields a reduction of polarization down to a ~ 0.5; still, 
this value corresponds to a significant transverse polar- 
ization. 

The asymmetry coefficient a was also computed for 
the directly produced J/ip and for the production via 
the xi decay [13]. The results are ot.<n r — 0.25 for di- 
rect production and a Xl ~ —0.15 for the production via 
Xi decay (except the forward region of xf > 0.8, where 
both otdir and a Xl begin to increase). After summing 
all channels of J/ip production it was found [13] that 
c? t e J t l ~ 0.5. Experimentally [17], no sizable polarization 
in the entire range of xp was observed, a ~ (there is 
however an indication that at very large xp a becomes 
negative). This disagreement between theory and experi- 
ment demonstrates again that the production mechanism 
of J/ip and possibly xi an d Xi m hadronic collisions is 
essentially non-perturbative. (Even though we have dis- 
cussed only ttN data, there is no reason to believe that 
in pN collisions the situation will be very different, apart 
from a relatively smaller contribution of the qq annihila- 
tion in the latter case.) It is also interesting to note that 
for the case of T production, the data from E866 Collab- 
oration [18] show transverse polarization for T(2S + 3S, 
in qualitative agreement with the predictions of pertur- 
bation theory. This of course is to be expected if the va- 
lidity of perturbation theory were to improve between the 
scales fixed by the masses of charm and bottom quarks. 

Let us now dwell upon the J/ip production in heavy 
ion collisions. Let us assume that at sufficiently high col- 
lision energy the quark-gluon plasma is formed. Due to 
the arguments presented above, the formation of quarko- 
nia will thus take place in the plasma. This will of course 
result in the suppression of the formation probability [1]; 
moreover, the presence of the plasma is likely to affect 
the excited states more significantly, and the contribu- 
tion of the excited quarkonium states to the observed 
yield of J/ip will thus change, which also can result in 
the change of the J/ip polarization [19]. If quarkonium 
is produced in the plasma, the non-perturbative effects 
should be absent (or small), and we are left only with 
the perturbative mechanism. Then, according to [13] 
about one half of J/ip : s will be produced directly and 
another one half via X2 — > J/ip + 7- (The approximate 
equality of these contributions stems from the fact that 
the extra power of a s in the direct production cross sec- 
tion is compensated by a relatively small branching ratio 



- about 20% - of the \2 —> J/ip + 7 decay; note also 
that x/J/ip ratio has been found to be independent of 
the collision energy - see [12].) We thus expect that the 
asymmetry coefficient of the electron (muon) angular dis- 
tribution in the J/ip — * e + e~(fi + fi~) decay in the case 
of quark-gluon plasma formation will increase from zero 
to about (at pt — 0) a ~ 0.6. The account of the initial 
transverse momentum distribution of gluons as discussed 
above reduces asymmetry coefficient to 



a ~ 0.35 -j- 0.4. 



(3) 



Still, we expect a remarkable increase in the asymmetry 
coefficient when going from hadron to heavy ion colli- 
sions. 

Of course, there are effects which may result in some 
decrease of a in comparison with (3), notably a more ac- 
curate account of the transverse momentum distributions 
of gluons and, as also discussed above, the interactions 
of J/ip with the constituents of hadronic and/or quark- 
gluon fireball. However, we do expect an increase of J / ip 
polarization in heavy ion collisions if the quark-gluon 
plasma is formed there. 
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